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ABSTRACT 

Our tightest upper limit on the sum of neutrino mass eigenvalues M v comes from 
cosmological observations that will improve substantially in the near future, enabling 
a detection. The combination of the Baryon Acoustic Oscillation feature measured 
from the Dark Energy Spectroscopic Instrument and a Stage-IV Cosmic Microwave 
Background experiment has been forecasted to achieve a(M„) < 1/3 of the lower 
limit on M v from atmospheric and solar neutrino oscillations (Abazajian et al. 2013; 
Fogli et al. 2012). Here we examine in detail the physical effects of neutrino mass on 
cosmological observables that make these constraints possible. We also consider how 
these constraints would be improved to ensure at least a 5cr detection. 

Key words: cosmology - cosmology:cosmic microwave background - cosmology: 
observations - large-scale structure of universe 


1 INTRODUCTION 

Basic questions about the neutrino mass matrix remain 
unanswered, such as whether the CP-violating phase is non¬ 
zero, whether the neutrinos are Majorana or Dirac, and 
whether the hierarchy of masses is normal or inverted (Fogli 
et al. 2012; Forero, Tortola & Valle 2012; Barger, Marfatia 
& Whisnant 2002). Significant experimental and observa¬ 
tional efforts are underway and being planned to answer 
these questions. Doing so may shed light on possible exten¬ 
sions beyond the standard model of particle physics. 

The question of the type of mass hierarchy may be set¬ 
tled by cosmological observations. The best lower limit on 
the sum of the mass eigenstate masses, M v , comes from 
analysis of solar and atmospheric neutrino oscillation data 
(Forero, Tortola & Valle 2014). If cosmological determina¬ 
tions of this quantity tighten up near this lower bound, then 
the inverted hierarchy will be ruled out. 

Current data lead to 0.058 eV < M„ < 0.21 eV, where 
the upper limit comes from cosmic microwave background 
(CMB) and baryon acoustic oscillation (BAO) data (Planck 
Collaboration XIII 2015). Bringing this upper limit down is 
a major science goal of a Stage-IV (S4) cosmic microwave 
background project, CMB-S4, and also of the galaxy sur¬ 
vey project Dark Energy Spectroscopic Instrument (DESI). 
These projects are forecasted to determine M v with a one- 
standard deviation of 45 meV (CMB-S4 alone) and 16 meV 
(CMB-S4 combined with DESI BAO) (Abazajian et al. 
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2013). These uncertainties are small enough to guarantee 
a detection of / 0 at a 3a or greater confidence level. 

Of course any conclusions from such cosmological data 
will be model dependent. How convincing will these data 
be that we are indeed seeing the impact of neutrino mass, 
and not misinterpreting some other signals? The forecasted 
precision is also not quite as strong as one would like; is there 
a way to guarantee at least a 5cr detection of M„ ^ 0? Here 
we address these questions. To address the first question 
we examine the particular signatures of neutrino mass that 
lead to the above forecasts. To address the second we look 
at what additional types of data can further tighten the 
expected uncertainties. 

The paper is organized as follows. In Section 2, we 
briefly introduce the cosmological signatures of massive neu¬ 
trinos. In Section 3, we focus on changes in the cosmic expan¬ 
sion rate and structure growth rate due to massive neutrinos. 
In Section 4, we analyze the influence of massive neutrinos 
on the CMB lensing potential power spectrum. Forecasts on 
the constraints of total neutrino mass from CMB and Large 
Scale Structure (LSS) measurements are given in Section 5 
and conclusions are presented in Section 6. 


2 SIGNATURES OF NON-ZERO NEUTRINO 
MASS 

The cosmological signatures of massive neutrino have been 
investigated since decades ago, e.g. (Doroshkevich et al. 
1980; Doroshkevich & Khlopov 1981; Doroshkevich et al. 
1981). We can more broadly view the cosmological neutrino 
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program as a study of the dark radiation that we know exists 
as a thermal relic of the big bang. By dark radiation here we 
mean anything, other than photons, thermally produced in 
the early universe that is relativistic at least through decou¬ 
pling. We know that such a background of nearly massless 
non-photon radiation exists with high confidence from light 
element abundances and the cosmic microwave background 
damping tail. Both are sensitive to the history of the expan¬ 
sion rate, which depends on the mean density via the Fried¬ 
man equation. Combining Helium abundance and CMB data 
constrains the effective number of relativistic species to be 
lV e ff = 2.99 ± 0.39 (Planck Collaboration XIII 2015). 

Is this background entirely that of the 3 active neu¬ 
trino species? Is any part of it from something else? Could 
there be a significant excess of neutrinos over anti-neutrinos? 
These are interesting questions, also to be addressed by fu¬ 
ture CMB observations that will significantly tighten up con¬ 
straints on N e f{. Our confidence that the dark radiation is 
indeed that of cosmological neutrinos with phase-space dis¬ 
tributions as expected from the standard thermal history 
will be greatly increased if the constraints on N e s tighten 
up to a(N e ff) = 0.02, as forecasted, consistent with the ex¬ 
pected value of 3.046. For the purposes of this paper, we will 
assume this is what will happen. 

If we assume that the dark radiation background is en¬ 
tirely the active neutrinos with the expected phase space dis¬ 
tributions, the assumption of non-zero neutrino mass leads 
to very specific predictions for cosmological observables. 
First we consider the expansion rate as a function of red- 
shift. The rest-mass energy of the neutrinos begins to slow 
down the decline of energy density with expansion as they 
become non-relativistic, leading to an increase in H(z) rel¬ 
ative to the = 0 expectation. This increase would per¬ 
sist to a = 0 if we were holding the other contents of the 
low-redshift universe constant. However, for our purposes of 
exploring observable consequences of M v ^ 0 it makes much 
more sense to hold the angular size of the sound horizon on 
the last-scattering surface constant, since this quantity is so 
well-determined from CMB observations (Planck Collabo¬ 
ration XIII 2015; Hinshaw et al. 2013). To do so one must 
decrease the density of dark energy. Assuming the dark en¬ 
ergy is a cosmological constant, the shape of A H(z) has a 
very particular form, as shown in Fig. 1, changing sign very 
near 2 = 1 , with the onset of dark energy domination. 

Were we able to trace out this departure of H ( 2 ) from 
the M„ = 0 shape, it would contribute to our confidence 
we are seeing the impact of non-zero neutrino mass. How¬ 
ever, as we will see, the DESI determinations of H(z) will 
be insufficient to resolve this very small signal across red- 
shift. That is not to say the signal is altogether observably 
invisible. These changes to H ( 2 ) affect comoving angular di¬ 
ameter distance Da(z) in ways that are detectable by DESI. 
It is just that it will be difficult, if not impossible, to make 
the case that there is the sign change in the H ( 2 ) correction 
near 2 ~ 1. The changes to H(z) also directly impact the 
growth of structure, with observable consequences for the 
redshift-space distortion (RSD) and the CMB lensing po¬ 
tential power spectrum, which we will discuss in Section 3 
and 4 respectively. 





Figure 1. The dependence of expansion rate H(z) and comoving 
angular diameter distance Da(z) on M u , where we minimize the 
X 2 (©, M„) by adjust the 6 ACDM parameters 0 when increasing 
M„ from 0 to 50,100, 200 meV. 


3 INFLUENCE OF MASSIVE NEUTRINOS ON 
GALAXY SURVEY OBSERVABLES 

To quantify the influence of massive neutrinos, we compare 
a fiducial cosmology with massless neutrinos and a cosmol¬ 
ogy with massive neutrinos. The fiducial cosmology is a 
flat ACDM universe with the Planck best fit parameters 
(Planck Collaboration XVI 2014), i.e. uJb = 0.022032, uj m = 
0.14305, A s = 2.215 x 10" 9 ,n s = 0.9619, r = 0.0925, H 0 = 
67.04 km/s/Mpc, M„ = 0 meV. The set of M„ 7 ^ 0 cos¬ 
mologies have parameters 9 = (0, M„), where 0 are ACDM 
parameters. Given a specific M v , we choose © by minimizing 
X 2 (0,M„), with 

X{Q,M V ) = FapXaX/3 

= FijXiXj + 2 Fi „ A t M „ -f- PM v , (1) 

where F is the Fisher matrix for the CMB observations, 
Xi = (0 — 0fld)i with i indexing the 6 ACDM parameters 
and summation over repeated indexes a,/3,i,j is implied. 
Minimizing y 2 (0,M l ,) requires 

0 = 8 X 2 /dXi -0 Xi = (2) 

where G is a subset of the Fisher matrix F, Gij = Fij. 

In Fig. 1, we show the influence of M u = 50,100 and 200 
meV on expansion rate H(z) and comoving angular diameter 
distance Da{z). We see that H(z < 1) decreases and H(z > 
1 ) increases compared to the fiducial cosmology, and the 
comving distance Da{z) increases accordingly. Though the 
departure of H(z) from the M v = 0 shape is undetectable 
by DESI BAO, the changes in Da{z) are readily detectable 
(Font-Ribera et al. 2014). 

One of the observable consequences of these changes to 
H(z) is the impact on the structure growth rate. How one 
describes this impact depends on what one is using for a 
comparison model. We use, as a comparison model, a cos¬ 
mology with massless neutrinos in place of the massive ones. 
One could also use as a comparison model one with addi¬ 
tional cold dark matter in place of the neutrinos. We use 
the former, consistent with our underlying assumption that 


© 0000 RAS, MNRAS 000, 000-000 


























































Constraints on neutrino mass from Cosmic Microwave Background and Large Scale Structure 3 



Figure 2. The dependence of structure growth rate on M„, where 
we minimize the x 2 ((-). AU) by adjust the 6 ACDM parameters 
© when increasing M„ from 0 to 20 meV. 


sition of structure growth rate from suppressed to enhanced 
is delayed to z ~ 0.5 compared to z — 1 when the expansion 
rate transitions. The reason for this delay is that the growth 
rate dD e (a)/da is determined by two factors: the expan¬ 
sion rate and gravitational attraction. The slower growth 
at z < 1 leads to weaker gravitational potentials. Around 
z ~ 1, the expansion rate is the same for the two models, 
but the gravitational potential is weaker for the M„ 0 
model. Therefore the growth rate remains suppressed,until 
some later time z ~ 0.5 when the weaker gravitational po¬ 
tential is compensated by even slower expansion. 

DESI will provide a comprehensive survey of spectro¬ 
scopic galaxies and quasars covering redshifts 0.1 < z < 3.5, 
with precision in each redshift bin Az = 0.1 better than 
a f<r S (a)/ fa's(a) = 0.35% from 0.4 < z < 1.5 (Huterer et al. 
2013). In Fig. 2, we show the influence of M v = 20 meV on 
the structure growth rate f(a)ag(a), which is detectable by 
DESI RSD. 


we have 3 neutrino species with phase-space distributions as 
expected from the standard thermal history. 

The structure growth rate that is usually quantified by 
d In erg (a)/d In a, can be determined by RSD from galaxy sur¬ 
veys, where crg(a) is the amplitude of mass fluctuations a a 
on scale of 8/V 1 Mpc, i.e., ag(a) = a(R = 8/i -1 Mpc;a). 
Here 

/ °° 7,3 

^Ps{k-,a)W 2 (kR)d\nk, (3) 

where Ps(k,a) is the matter power spectrum defined by, 
(5(k;a)<5*(k';a)> = P s (fc; a)d (3) (k - k'), and W(kR) = 
3ji(kR)/kR. is the window function. Introducing the growth 
function D(a), 


D(a) 


S(a) 

S(a = 1) 


o~8 (a) 

ag(a = 1) ’ 


( 4 ) 


where 5{a) is the matter overdensity at redshift z = 
1/a — 1, and defining perturbation growth rate f(a) = 
d In D(a )/d In a. (Pierpaoli, Scott & White 2001), we may 
rewrite the structure growth rate as d In <78 (a)/d In a = 
f(a)a 8 (a). 

In the above definition, the growth function D(a) and 
the structure growth rate f(a)ag(a) depend on ag(a = 1) 
which varies for our fiducial M„ = 0 model and the model 
with massive neutrinos. For easier comparison, we introduce 
the growth function D e (a) normalized at early time, say at 
a e = 1/1100 when massive neutrinos are relativistic, 


D e (a) 


^( a ) 

5(a e ) 


os (a) 
ag (a e )’ 


( 5 ) 


and the structure growth rate can be rewritten as 

/(a)ag(a) = a s(a e ), (6) 


where ag(a e ) is the same for the model with massive neu¬ 
trinos and the fiducial model. So the variation of f(a)ag(a) 
only depends on the growth rate dD e (a)/da. 

In Fig. 2, we show the impact of massive neutrinos with 
mass 20 meV on the structure growth rate. The structure 
growth rate is decreased at high redshift because of the en¬ 
hanced expansion rate (Fig. 1). Note though that the tran- 


4 INFLUENCE OF MASSIVE NEUTRINOS ON 
THE CMB LENSING POWER SPECTRUM 


4.1 Introduction to the lensing power spectrum 


We begin with a brief review of gravitational lensing of the 
CMB. For details see e.g. the review by Lewis & Challinor 
(2006). Gradients in the gravitational potential, 4’, distort 
the trajectories of photons traveling to us from the last scat¬ 
tering surface. The deflection angles, in Born approximation, 
are d = V0, where the lensing potential, cf>, is a weighted 
radial projection of <f>. The key quantity for calculating the 
impact of lensing on the temperature power spectrum is the 
angular power spectrum of the projected potential, Cf^, 
which we also call the lensing power spectrum. Taking ad¬ 
vantage of the Limber approximation (Limber 1953), it can 
be written as a radial integral over the three dimensional 
gravitational potential power spectrum P<i> 


1 4 C?* ~ 4^*‘ d X {k 4 Pf) 



( 7 ) 


where x is the comoving distance from the observer, a = 
a(x), a * subscript indicates the last scattering surface, (1 — 
x/x*) 2 > s the lensing kernel, and the power spectrum is 
defined as 


($(k; a)4»*(k'; a)) = P<s>(fc; a)<S (3) (k - k'). (8) 

To calculate Pi. we assume a power-law primordial power 
spectrum P|(fc), 

k 3 { k \ n » _1 

^PJW^(A) , ( 9 ) 

where ko is an arbitrary pivot point, A s and n s are the 
primordial amplitude and power law index respectively. 

The gravitational potential at late times, <3?(A;, a), is re¬ 
lated to the primordial potential & p (k) by (e.g. Kodama & 
Sasaki 1984) 

4>(fc, a) = ^4> p (k)T(k)s(k ; a)g(a), (10) 

where the potential on very large scales is suppressed by a 
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Figure 3. The lensing power spectrum calculated from CLASS 
(solid line), calculated with Limber approximation (dashed line) 
and calculated with Limber approximation and setting g(a) = 1 
(dashed-dotted line). 

factor 9/10 through the transition from radiation domina¬ 
tion to matter domination. For modes which enter the hori¬ 
zon during radiation domination when the dominant com¬ 
ponent has significant pressure (p fts p/3) the amplitude of 
perturbations cannot grow and the expansion of the Uni¬ 
verse forces the potentials to decay. In a cosmology with 
massless neutrinos, for modes which enter the horizon after 
matter-radiation equality (but before dark energy domina¬ 
tion) the potentials remain constant. The transfer function 
T(k) takes this into account, being unity for very large scale 
modes and falling approximately as k~ 2 for small scales. The 
transfer function T(k) is independent of scale factor a be¬ 
cause in a cosmology with massless neutrinos, potentials on 
all scales keep constant in matter domination. The impact 
of massive neutrinos on the potentials is described by the 
function, s(fc; a) , which is unity on scales above the free- 
streaming scale and decreases with time on scales below. 
Once the cosmological constant starts to become important 
the potentials on all scales begin to decay. This effect is cap¬ 
tured by the growth function, g(a), which is unity during 
matter domination. With these definitions we have 

(fc 4 P$)(fc; a) = j^7rV(a) s 2 (fc;a) kT 2 (k) A s 

( 11 ) 

With these pieces in place, we now examine the accu¬ 
racy of the Limber approximation. According to Loverde 
& Afshordi (2008), it is a better approximation to replace 
k = £/x with y/e(t+ l)/x — (-£ + 0.5 )/x in the original Lim¬ 
ber approximation Eq.(7). With the replacement and defin¬ 
ing x = x/x *i the lensing power spectrum can be written 
as 

+ 1) 2 C^ ~ 4 X * [ dx (fc 4 P<i>) A (i _ 

Jo \ X X* J 

( 12 ) 

where a = a(x). Fig. 3 shows the lensing power spectrum 
calculated from CLASS (Lesgourgues 201 la, b; Bias, Lesgour- 
gues & Tram 2011; Lesgourgues & Tram 2011) compared to 
the lensing power spectrum calculated from the Limber ap¬ 
proximation. We see the Limber approximation reproduces 


Figure 4. The dependence of the lensing power spectrum on total 
neutrino mass, (Cf^ — C/)/ ] )/C| , / (1 = Re X (M„/eV). The black 
line is the numerical result from CLASS, the red line is the result of 
Limber approximation setting x* = 1.4 X 10 4 Mpc, the blue line 
is the result of Limber approximation with g(a) = 1, the magenta 
line is the result of Limber approximation fixing both ui m and 
g (a), and green line is the result of Limber approximation with 
A 3 ,n 3 , a; m ,tuand g(a) fixed. 

the accurate numerical result even for small £. In order to 
understand the the influence of the growth function, we also 
calculated the lensing power spectrum by setting g(a) = 1. 
The growth function makes a difference for £ < 50 (Pan, 
Knox & White 2014). Now we are to understand the impact 
of massive neutrinos on Cf^. 

4.2 Influence of massive neutrinos on the lensing 
power spectrum: results 

According to Eq.(ll) and Eq.(12), it is clear that Cf^ is de¬ 
termined by the primordial perturbation A s , n B , the transfer 
function T(k), the impact of massive neutrinos s(k;a), the 
growth factor g(a), and the comoving distance to the last 
scattering surface 

In order to quantify the dependence of the lensing power 
on total neutrino mass M„, we take samples from a Planck 
ACDM +M„ chain, and fit the linear relation 

('•«£'#’ _ 

-—^-= R e x M v (eV), (13) 

( ~'e,Rd 

The linear fitting result is shown as the thin solid line in Fig. 
4. To understand the contribution of the various parameter 
variations and effects, we also plot Re in Fig. 4 for the cases 
shown. 

We now work our way towards an understanding of the 
full response, the thin solid curve, in stages. We begin with 
the case where we fix most of the parameters other than 
M„, and turn off the impact of dark energy on the growth 
factor, g(a), fixing it to unity. In this case, increasing M„ 
has two effects: 1) it decreases the free-streaming length 
~ ( T v /Mu)/H{z ) (Abazajian et al. 2011) , and 2) increases 
the expansion rate once the neutrinos start to become non- 
relativistic. The increased expansion rate acts to suppress 
the growth of structure on all length scales. However, the de- 
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creased free-streaming length acts to boost structure growth 
on scales above the free-streaming length, nearly exactly 
canceling the suppression. The result is the bottom-most 
curve of Fig. 4: nearly no effect at low £, a constant suppres¬ 
sion of power at high l, and a smooth transition between 
these two regimes. 

The difference between the bottom-most curve and the 
dot-dashed curve (the one labeled, ‘Limber with g(a) = 1’) 
is due to how other parameters adjust as M v varies. We can 
isolate these changes as almost all due (at least at £ > 50) 
to a correlation between M„ and u) m , as demonstrated by 
the following: if we fix cn m , and let A s , n s , and Ub vary, 
we get the second curve from the bottom which differs very 
little from the bottom-most curve. Once we let uim vary as 
well, we get the dot-dashed curve. Letting vary leads to 
variation in n s and A s as well, but these changes all flow 
from the correlation between uj rn and M„. 

Once we let g(a) vary as well we get a curve that is in¬ 
distinguishable from the thin solid curve, which is boosted 
everywhere, and especially at £ < 50. The contributions to 
these large angular scales come predominantly from modes 
to the left of the peak in the matter power spectrum. At 
fixed (large) angular scale, structures that are nearer by, and 
therefore on smaller length scales, are closer to the peak of 
the matter power spectrum. Thus the large angular scales 
are weighted toward later times, and therefore more influ¬ 
enced by g(a ) than the smaller angular scales. The growth 
factor increases with increasing because to keep the an¬ 
gular size of the sound horizon fixed Cl a must decrease. 

The lensing kernel’s dependence on cosmological param¬ 
eters comes entirely via its dependence on x*- To see how 
much of the variation in the lensing power spectrum is due 
to the lensing kernel, we fix x* = 1-4 x 10 4 Mpc (top-most 
dashed curve). By examining the difference between the top¬ 
most dashed curve, for which x* is fixed, and the thin solid 
curve, which is the full numerical result, one can see this 
effect is very small 1 


1 This result is in contrast to the case of tomographic cosmic 
shear as a probe of dark energy. In this case the sensitivity of the 
data to variations in the dark energy equation-of-state parameter 
largely arises from the lensing kernel (Simpson & Bridle 2005; 
Zhang, Hui &; Stebbins 2005). 

To summarize, there are three main effects of massive neutri¬ 
nos on the lensing power: 1) increased expansion rate suppresses 
power, 2) decreased free-streaming length compensates for the 
suppressed power at scales above the free-streaming length, 3) 
other parameter variations due to partial degeneracies in C'f 1 
(most notably an increase in u} m ) boost the power on all scales. 
The net result is increased power at large scales and a decrease in 
power at small scales. One might potentially include the growth 
factor here as the fourth-most important effect, somewhat in¬ 
creasing the power at large angular scales. 

The origin of the degeneracy in Cf T between and M v is 
actually due to lensing itself. Planck Collaboration XVI (2014) 
demonstrated that the dominant effect leading to the constraint 
of neutrino mass from the CMB temperature anisotropy power 
spectrum is gravitational lensing. As shown in Fig. 4, increasing 
suppresses the lensing power, while increasing iu m increases 
the lensing power. The lensing power suppression by massive neu¬ 
trinos can be compensated by the enhancement from increasing 
<u m , so uncertainties in M„ and co m are expected to be positively 
correlated (Namikawa, Saito &c Taruya 2010). 


5 FORECAST OF CONSTRAINTS ON THE 
TOTAL NEUTRINO MASS FROM 
DIFFERENT DATA SETS 

We use the Fisher matrix formalism to forecast constraints 
on neutrino mass from future CMB and LSS experiments. 
The fiducial cosmology used here is the same as the one used 
in Section 3 except with a different value of total neutrino 
mass, M v = 85 meV. 


5.1 CMB-S4 and DESI BAO 


Following Wu et al. (2014) and Dodelson (2003), the Fisher 
matrix for cosmological parameters constrained by CMB 
spectra is written as 


F a p 


E 


i 


2^+1 

2 


/skyTr 


,-idCt_ r -idCA 
’ l de a l dOfi)’ 


(14) 


and it is related to the expected uncertainty of a parameter 
9 a by a(d a ) = ■\/(F’- 1 ) aQ ,, where 

( cJ T + nJ t cJ e cj d \ 

C e = Cj E Cf E + Nf E 0 

V Cj d 0 Ci d + Nf d ) 

(15) 

and C dd is the angular power spectrum of the deflection 
field d, which is related to the lensing power spectrum by 
C dd = £(£ + 1)C/^. The Gaussian noise N xx is defined as 

N t xx = A 2 x exp (l(£ + 1) ) , (16) 

where A.y (X = T, E, B ) is the pixel noise level of the ex¬ 
periment and #fwhm is the full-width-half-maximum beam 
size in radians (Knox 1995; Zaldarriaga, Spergel & Seljak 
1997). The noise power spectrum of deflection field Nf d is 
calculated assuming a lensing reconstruction that uses the 
quadratic EB estimator(Okamoto & Hu 2003). We use the 
iterative method proposed by Smith et al. (2012), which 
performs significantly better than the uniterated quadratic 
estimators (Hirata & Seljak 2003). 

For the CMB-S4 experiment, we assume the temper¬ 
ature noise level At = 1.5 /rK-arcmin, the polarization 
noise level Ae = As = \/2 At, the fraction of covered 
sky / a ky = 0.5 and the beam size #fwhm = l r . With these 
given experiment sensitivities, we obtain a constraint from 
CMB with a(Mv) = 38 meV. The 1 <t and 2<r constraint are 
shown in Fig. 5. 

According to the analysis in Section 3 and 4 , DESI 
BAO are helpful to break the degeneracy between Mv and 
cu m . BAO uncertainties are independent from CMB experi¬ 
ments, so the total Fisher matrix is simply given by addition 

Fcmb+bao = Fcmb + Fbao, (17) 

where the DESI sensitivities of BAO signal can be found in 
Font-Ribera et al. (2014) and shown in Fig. 1. It is found 
that, adding the DESI BAO data greatly improves the con¬ 
straint to a(Mv) = 15 meV (similar forecasts were also con¬ 
ducted by Abazajian et al. (2013); Wu et al. (2014)). 


5.2 Beyond DESI BAO 

The largest signal of massive neutrinos on H(z) and D(z ) is 
found at low redshifts (see Fig. 1), where BAO has inevitably 
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M v [eV] 




Figure 5. Forecasted fcr and 2 a constraints in the AC — oj m 
plane, where the CMB-S4 experiment results in a ct(M„) = 38 
meV constraint, the combination of CMB-S4 and DESI BAO 
yield a <r(M„) = 15 meV constraint, and adding measurements 
of the structure growth rate by DESI RSD further improves the 
constraint to cr(M„) = 9 meV. 

large noise because of small amount of survey volume and 
large cosmic variance. Other than DESI BAO, we also inves¬ 
tigate other low-redshift tracers of H(z) and D(z) which are 
possible to tighten the uncertainty of total neutrino mass. 

DESI RSD: similar to BAO, RSD uncertainties are also 
independent from those of CMB observations, so the total 
Fisher matrix of CMB+BAO+RSD is also approximately 
given by addition 

Ecmb+bao+rsd = Ecmb + Fbao + Ersd, (18) 

where we use the RSD sensitivities from DESI survey which 
can be found in Huterer et al. (2013) and shown in Fig. 2. 
Here we use the approximation that uncertainties in BAO 
and RSD are uncorrelated, due to they are sensitive to differ¬ 
ent aspects of the matter power spectrum: BAO is sensitive 
to its characteristic length scale r s while RSD is sensitive 
to its amplitude. In fact, our result is insensitive to the 
approximation because we find that both CMB-S4+DESI 
BAO+DESI RSD and CMB-S4+DESI RSD yield the same 
cr(M„) = 9 meV uncertainty. 

Better BAO: DESI survey cover 14, 000 squared degrees 
(about 1/3 of the whole sky). We explore a future BAO 
experiment which covers the whole sky and in which cosmic 
variance dominates over shot noise in the redshift range 0 < 
« < 4.0. Constraints on Da{z) and H(z ) from this BAO 
experiment are shown in Fig. 6. It is found that CMB-S4 
and the cosmic variance limited BAO constrain the total 
neutrino mass with uncertainty a(M v ) = 11 meV. So we 
conclude that 11 meV is a lower limit of a{M v ) we could 
measure from CMB-S4+BAO, where the limit mainly comes 
from noise level of the CMB lensing signal. 

Supernovae: The constraining power of BAO is limited 
by its large cosmic variance at low redshifts (Fig. 6), so su¬ 
pernovae distance measurements which do not suffer from 
the cosmic variance problem may be effective complements 
if their systematic errors are well controlled. Supernovae per¬ 
form better in relative distance measurements than in abso¬ 
lute distance measurements. However for the ACDM + M„ 


Figure 6. Same as Fig. 1, but with suppressed errorbars of D.\ (z) 
and H{z) coming from CMB-S4 and a cosmic-variance-limited 
BAO experiment. 



Figure 7. The uncertainties in relative distances from CMB-S4 + 
DESI BAO. Note that the uncertainties is multiplied by a factor 
of 10 4 in the plot. 

model, the uncertainties in relative distances from CMB-S4 
+ DESI BAO are very small (see Fig. 7) . We conclude that 
supernova observations must result in relative distance de¬ 
terminations with systematic errors less than about 0.05% if 
they are to tighten the constraints on neutrino mass. Com¬ 
pared to systematic errors from current supernova observa¬ 
tions (e.g., Suzuki et al. (2012)) this would be a reduction 
by a factor of ~ 20 . 


6 CONCLUSION 

This paper is motivated by our desire to better understand 
the origin of current and forecasted cosmological constraints 
on the sum of neutrino masses. We took as a given that de¬ 
termination of N e g will solidify the predicted value of 3.046, 
increasing our confidence that the phase-space distribution 
of the cosmic neutrino background is what we expect based 
on the standard thermal history. With that as a given, the 
most important aspect of increased neutrino mass (relative 
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to some reference model) is an increased neutrino energy 
density. If the model with increased mass is to remain consis¬ 
tent with CMB observations, the distance to last-scattering 
must be preserved and so the total energy density, and there¬ 
fore the expansion rate, cannot increase at all redshifts. To 
compensate for the increase in neutrino energy density, the 
cosmological constant must decrease in value. Thus varying 
neutrino mass leads to changes in H (z) with a very partic¬ 
ular shape: a mild increase at high redshifts, a larger de¬ 
crease in low redshifts, with a transition near the onset of 
A—domination at 2 ~ 1. Unfortunately this very specific 
prediction of the shape for H(z) is difficult to verify in de¬ 
tail because of how small the departures from the reference 
model are at z > 1. We see that this difficulty persists even 
for a cosmic-variance limited all-sky (z < 4) BAO experi¬ 
ment (see Fig. 6). 

The sensitivity of CMB-S4 to neutrino mass comes via 
the impact of this increased expansion rate on the growth 
of structure. At scales below the neutrino free-streaming 
length, this increased expansion rate suppresses the growth 
of structure. Above the free-streaming length, the ability of 
massive neutrinos to cluster compensates for the increased 
expansion and there is no net suppression. Because increas¬ 
ing matter density increases lensing power amplitude, the 
CMB lensing-derived constraints on neutrino mass have un¬ 
certainties positively correlated with the matter density un¬ 
certainties. This correlation with dark matter density leads 
to secondary correlations of neutrino mass uncertainty with 
uncertainties in n s and A B . We disentangled all these various 
effects in Fig. 4. 

The correlation between and uj m has the opposite 
sign as that from BAO, since increasing neutrino mass and 
increasing ui m both increase the expansion rate at z > 1, and 
lead to a compensating decrease at 2 < 1. Thus the com¬ 
bination of CMB-S4 and DESI-BAO leads to improvements 
in the determination of both quantities. 

Finally, we briefly investigated how constraints might 
be improved beyond the ~ 3a to ~ 4a detection expected 
from CMB-S4 + DESI BAO in the case of the lowest possi¬ 
ble neutrino mass of 58 meV. The larger signals are at low 
redshift, so we considered supernovae as relative distance 
indicators and redshift-space distortions, as well as cosmic- 
variance-limited BAO. We found requirements on supernova 
precision that are probably prohibitively stringent. But bet¬ 
ter BAO, and RSD, both have the potential to improve the 
detection to ~ 5a or greater. 
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